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INTRODUCTION 


The interaction of elastic waves with material inhomogeneities 
such as second -phase particles and cracks plays an important role in 
material and fracture characterization by quantitative ultrasonics [1,2]. 
These inhomogeneities represent mis-matches in material moduli as well 
as mass density compared with the surrounding matrix material and thereby 
generate elastodynamic fields in addition to the fields when they are 
absent. These fields are referred to as the scattered fields when mea- 
sured far from the inhomogeneities. When the non-dimensional wavenumber 
is greater than unity, ka > 1 , the term diffraction is usually used in 
place of scattering. 

Diffraction of elastic waves by thin flat objects, cracks in parti- 
cular, has been an actively pursued subject for sometime by theoretical 
mechanicians. Due to the complexity of the presence of a third dimension, 
most work is confined to essentially a two-dimensional situation. Re- 
cently, Teitel [3], Gubernatis and Domany [4], and Datta [5], using the 
quasi-static approach and asymptotic expansion approach, respectively, 
obtained results for an elliptical crack that are valid in the Rayleigh 
limit, ka<<l. Budiansky and Rice [6] gave a general integral formulation 
for the dynamic response of an isolated three-dimensional crack. Employing 
the geometric diffraction theory [7], Gautesan, Achenbach and McMaken [8] 
investigated the diffraction by elliptical cracks valid for the region of 
short wavelength, ka>>l. A summary of the comparison with experimental 
results is given in [9] for a penny- shaped crack and [10] for an elliptical 
crack. For the region of medium wavelength, results available are obtained 
mostly by numerical methods through a matrix approach. 
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In the present paper the scattering of elastic waves by an isolated, 
flat, thin, and elliptical inhomogeneity is studied by employing the ex- 
tended version of the method of equivalent inclusion. The method of 
equivalent inclusion was originally developed by Eshelby [11,12] for 
the determination of elastostatic fields inside and outside an ellipsoidal 
inhomogeneity. Wheeler and Mura [13] first extended this method to study 
the elastodynamic response of composites where the mis-match in mass 
density is ignored. Recently, Fu and Mura [14] gave a complete formula- 
tion to extend the method of equivalent inclusion for elastodynamic problems. 
They presented the equivalence conditions that allow the elastodynamic 
fields inside and outside an inhomogeneity to be determined. In addition 
to the presence of eigenstrains , (1) , the concept of "eigenforces” ttJ 1 

is also introduced. 

The solution procedure given in [14] requires the evaluation of 
certain volume integrals that are associated with the inhomogeneous Helm- 
holtz equation. A method for evaluating these integrals is given by Fu 
and Mura [15]. It is easily seen from Ref. [15] that the results in Ref. [14] 
reduce to the elastoststic solution when the frequency approaches zero, 
co 0. It is further noted that since all the geometric information are 
contained in these volume integrals, scattering of inhomogeneity of any 
geometry can be obtained simply by evaluating these integrals. 

The case of the scattering of a perfect sphere is studied in detail 
by Sheu and Fu [16] and compared with classical results of Ying and Truell 
[17]. It was shown that the comparison is good up to ka about two when 
uniform distribution of eigenstrains and eigenforces are assumed. 


^Although this term was not used in Ref. [14] , it appears to be appropriate 
for this quantity as it can be seen from its definition = Apu... 

Its unit is of course the same as body force, i.e. force per unit volume. 
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In this paper the solution for a flat, thin elliptical inhomo- 
geneity is obtained by collapsing the ellipsoid, say letting a 3 -*■ 0 
Since the eigenstrains and eigenforces become infinite the limiting 
concept described by Mura [18] is employed. 
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ELASTODYNAMIC FIELDS OF AN ISOLATED INHOMOGENEITY 


For convenience in reading, it appears to be worthwhile to summarize 
the solution given in Fu and Mura [14], Consider the physical problem 
of an isolated inhomogeneity embedded in an infinite elastic solid which 
is subjected to a plane time-harmonic incident wave field as depicted 
in Fig. 1. Replacing the inhomogeneity with the same material as that of 
the surrounding medium, with moduli Cj^ rs and mass density p, and include 
in this region a distribution of eigenstrains and eigenforces, the physical 
problem is now replaced by the equivalent inclusion problem. 

Following Fu and Mura [14] the total field is now obtained as the 
superposition of the incident field and the field induced by the presence 
of the mis-matches in moduli and in mass density written in terms of 

*(13 


eigenstrains and eigenforces, nt: 


F = F^ + F^ 


( 1 ) 


where F denotes either the displacement field u^ , the strain field , 
or the stress field . The superscripts (i) and (m) denote "incident" 
and "mis-match", respectively. 

For uniform distributions of eigenstrains and eigenforces, the fields 
can be obtained as: 


uW(f, - 


( 1 ). 


- IT* S. (r) - C-, e* 'S . , 
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mn 
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where a comma denotes partial differentiation and 


S.(r) = / g. (?-*') dV 


jm 




( 5 ) 


in which g^ is the spatial part of the free space Green f s function and 
ft is the region occupied by the inhomogeneity. It is noted that the 
integrals S ^ and their derivatives must be evaluated for the region 
r > r' and for the region r < r’. Ref. [15]. The solution form repre- 
sented in Eqs. (1-5) gives the fields inside and outside an isolated 
inhomogeneity of arbitrary shape. 
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AN ISOLATED FLAT THIN INHOMOGENEITY 


Let the incident displacement field be longitudinal and of frequency 


a), amplitude u q : 

= U D qj expfiox^ - iut] (6; 

2 

where i = -1 and is the unit vector in the normal direction of the 
plane time harmonic wave. For a linear isotropic medium, the spatial part 
of the free space Green's function is well known. Substituting g^Cr-r') 
in Eq. (2) and using the limiting concept* 

limit a„ 7r* = A. constants (7) 

V° ' ' 

limit a_ = B. . constants (8) 

V° 3 « 


fs 1 — 

the scattered displacement lT j (r,t) from a thin elliptical flat inhomo- 
geneity can easily be obtained as: 


l“V 3 u o 


< oa l> 3u 0 r 


= [(C G m expiar)/ar + (D H m expiBr)/Br] exp(-io)t) 


G »- - <.^ a inVj^ c,MH1 - 2a /B )vjj* 2(a ' /6 ' ^VkVk 1 

H m ' ta 2 /. i nCE/a) 3 C^^ J -6 m;( 3AJ( 1 p/p) 

- 2 «'“>V{. * 2 ( 6 /“) 2 Wj B kj) 


l This concept is discussed in [18] in the light of Eq. (8). 
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c = 


oo 


I 

n=0,2 


n (-l) n C<ia 1 ) n (a 2 /a 1 ) n - A Aj^- ,! 'Cn/2) ! 
gl 0 '2 (n+3)(n+l)n! (fc/2) ! (n-fc)/2! 


D 


oo 

I 

n=0,2 


I 

J> = 0,2 


(-l) n (3a 1 ) n (a 2 /a 1 ) n *^^ _A Cn/2) ! 
(n^3T (n+l)n! (*/2) ! (n-*)/2! 


(12) 


(13) 


in which m, j ,k=l , 2 ,3, and (a^ > a 2 ), Z^, a, 3, and A p denote the semi-axes 

of elliptical inhomogeneity, direction cosines of scattered displacements, 

longitudinal wavenumber, shear wavenumber, and mass density mis-match 

(p’-p)* respectively. Also, At and Bt, are the reduced non-dimensional form 

3 3 K 

of Aj and B^ k , respectively, defined as follows: 

AJ - A^/dpoi 2 ^) (14) 

B )k ■ «« 

Expressions for the differential cross section dP(to)/dft and total 


cross section P(o>) can be obtained as before 

= a p (e ,<t>) + (a/ 6)a S (6 ) (16) 

P(u>) = / a p (e,4>) + (ci/B)ct s (0 ,<f>) df2 (17) 

Z 

where dft is the differential element of solid angle and 

a 2 a p (e,<j») = (ccaj) 6 (C GJ [C“G^] (18) 

e 2 a s (e,«j>) = ( aai ) 6 [D H m ]["FlT] (19) 


in which the super bars denote complex conjugate. It is noted that the 

constants At and Bt, must be evaluated from the equivalence conditions 
3 3 K 

given in Ref. [14] with the use of the limiting concepts in Eqs. (7,8) 
and of the integration method developed in Ref. [15] . 
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DETERMINATION OF A* AND B* 

I Jk 

In Eqs. (9-13) the scattered displacement field is given in terms 

of the "reduced" form of the eigenforces and eigenstrains, i.e. At and 

Bt, , see Eqs. (7,8,14,15). These constants must in turn be determined 
J K 

from the equivalence conditions, Eq. (14), Ref. [14]. Writing the incident 

wave field in a Taylor series the governing simultaneous algebraic equations 

can be easily obtained. Since f ..[0] and F ..[01 vanish automatically, 

these governing equations become uncoupled and lead to a three by three 

system for At and a six by six system for Bt, , where Eqs. (25,26) in 
J J K 

Ref. [14] are used. For a linear elastic medium, they are: 


Apoo u f. [0]At + A* = - q 
O JS l J j s n s 


( 20 ) 


{ AX6 D .. [0] + 2 ApD [0] }B* + { A6 B* + 2pB* 

st mmjk L J stjk L J jk st mm H st 


} 


= - + 2 Aw V*t } (21) 

where the subscripts s,t,m,j,k = 1,2,3, repeated subscripts denote sum 
from 1 to 3, and 


2 2 

4ttpo) f. (r) = - g <p 6. + iK • - (f>, . 

K js K J Y js mj 


( 22 ) 


4 " pu D stjktf> = 2p[ ^stjk + ♦’stjkl 


- fB [4>, jt 6 ks + ^ js <5 kt ] 


- Xa \p, 6 . , 

v mn jk 


(23) 


AX = X ' - X , Ay = F ' ~ F » Ap = p ' - p 
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I 


The <f>- and ^-integrals and their derivatives are evaluated by 
the method suggested in Ref. [ 15 ]. Retailing terms up to (aa.^ or (Ba^ 
of the fourth order, the constants are obtained as: 

2 

A? = - qj/{u Q a) Apf \ [ 0 ] + 1 } , no sum on j 

f.[ 0 ] = (f n [ 0 ], f 22 to], £ 33 [ 0 ]) 

{Bf} = [b..]' 1 ^.] i , j = 1 > 2 , 3 

{Bt.} = - q. q •/ [l+SC^itOj+ 4 *. [ 0 ])] , no sum on i,j . 

1J xjtj 3 13 

i,j=l,2,3 C = 1/4 tt 

where in Eq. ( 25 ) B* = BJ^ , B* = B 22 , Bj *= B*^ , 
c^. = - ( AA +2 Ap) [1 + qj] , j = l, 2 , 3 , 

b ll = ( A+2y ) + XCt.jjlO] + 2 (A+ 2 p) jj [ 0 ] + 2 yC'l', 11 [ 0 ] 

b 12 = A + 2 + 2£^<f> 22 [0] 

b 13 = A + + 2 pCV '»21 [ 0 ] + 2 X£;<}> , 22 [ 0 ] 

b 21 = A + + 2u£i{/, 22 [0] + 

b 22 = ^ A+2lJ ^ + H^jj[0] + 2\j£\p, 22 [0] + 2C(A+2u)4>, 22 [0] 

b 23 = A + + 2vjC'J j > 22 [0] + 2£A<J>, 33(C)] 

b 21 = A + A^jjtO] + 2 vt\l>, 33 [ 0 ] + 2 XQ<p, ^ j [ 0 ] 

b 32 = A + Xc^»jj[0] + 2iiCr> 33 [0] + 2 \£<p, 22 [ 0 ] 

b 33 = (A+2y) + A^^fO] + 2y^, 33 [0] + 2 (> +2y) £ 4 >, 33 [0] 

S = 1/4 tt C = ( 1 - 2 q 2 /B 2 )/ 4 ti 


(24) 


(25) 

(26) 


(27) 
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Note that In Eq. (26), <J>j = . <t> 2 ~ *22^°^ * * 3 = ^33^°^ ' 

The f- , 4>- and ip- functions are given as: 


4tt PW 2 f js [0] = - B 2 <j>[0]6 js + ^ js [0] - *, js [0] (28) 

B 2 4> [0] = 7ra ia2 B 2 {I o - ( (Ba^ 2 /16) ^ + i (4/3) B(a 2 +a 2 ) } (29) 

*, n [0] = -[(Tra^a 2 B 4 )/12) • Ij (30) 

4>, 22 [0] = -[(Tm^B 4 )/^] • I 2 (31) 

4>, 33 [0] = 0 (32) 


in which 


i = / 


dip 

aOiO 


= F(0,k) 


i/^d iJj 


.,2 


/ 7 

o (a^+i/;) A (^) °1 k 


f- { E(ejO _ k F (e,k) } 
a, ,2 ,2 


tydtp = 2_ [ F(9,k) _ E(Q ,k) j 
o (a 2 -^)A(i^) a l k 2 k 2 


/ 


doj 


> ,..2.2 sl/2 

o (1-k sm u) 


, E = / (l-k 2 sin 2 w) 1/2 d< 


(33) 

(34) 

(35) 


and as a^, B+tt/ 2, k 2 -> (1 -a 2 /a 2 ) and k’ 2 = (l-k 2 )^a 2 /a 2 , if a x > a,. If 
a.+0 and a 1 =a 2> we have I = it , ^ = 1 2 =Tr / 2a i • 
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EXAMPLES: SCATTERING OF A CRACK AT LONGWAVE LIMIT 


Consider now the scattering of an elliptical crack as an example. 
Since the crack is a void, 

Ap = -p , AA = -A , Au = -y 

The constants A? and B* can be easily obtained by solving the simul- 

j j k 

taneous equations, Eqs. (24-26): 

A j = Cl-ai a 2 ^ 2;p C e (36) 

B ij = _q i q j + + if i^ j > no sum on i,j. (37) 

B ij = -( 1+ qiq-jHU + 2u) (-X+2u)]/2v J (3X + 2 V i) 

+ ^i3^j3 [(* + 2y) (2A) ] /2p (3A+2y) if i=j, no (38) 

sum on i , j . 

For ultrasonic applications, the longwave limit is always of interest, 
i.e. only terms up to the second order in (aa) and (ga) are retained in 
the scattered displacement amplitudes. From Eqs. (9-13), it is observed 
that, this requires that only terms of (aa) and (ga) to the zeroth order 
be retained in Eqs. (36-38). A direct substitution of Eqs. (36-38) into 
Eqs. (9-13) easily yields the expression for the scattered displacement 
amplitudes in the longwave limit. Simple manipulation yields: 

CG m = (a 2 /3a 1 )Jlm{-£ jqj + (a 2 /6 2 )(£ j a kq .q k ) jA=lj2> 3 

- (2- 3 2 / a 2 ) (2-2v-vq^)/(l+v) (l-2v) 

+ ( 1 + ^q k )U-3v)/( 1 + v) +(l + q 2 )£ 2 (2v)/(l + v)} (39) 

and 
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DH m - (a 2 /3a 1 )((6/c,) 3 ( t]1 , tr4m;i )q j 
* (6/a) (-(*„, 

- (6/a) 4 Ji m [(l-(£ 1 + £ 2+ £ 3 )-Jl j q j 2 + £^) (l-3v)/ (1+v) 

+ (l + q 2 )()l 2 -£ 3 )(2v)/(l + v)]} (40) 

where v is Poisson’s ratio. 

Let the incident wave be going from the negative z-axis to the 
positive z-axis, Fig. 2, it is clear that 
qi=q 2 = 0 q 3 =l 

For the forward and backward scattered displacements along the z-axis, 
the direction cosines are £^=£ 2=0 and £^=+1, respectively. 

The displacement amplitude can be found, for example, as follows: 

u^ S ^ = u^ S ^ = 0 (41a) 

u^ S V (aaj)^^ = (a 2 /3a 1 ) { -1 +2 (1 -v- 3v 2 ) / (1 + v) (1 -2v) } (expiaz)/az (41b) 

where the upper and lower signs refer to the forward and backward scattering 
along the _+z-axis, respectively. 
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CLOSING REMARKS 


The solution to the direct scattering of an elliptical inhomogeneity 
is obtained by employing the extended version of Eshelby's method of 
equivalent inclusion and a limiting concept. When the inhomogeneity becomes 
a flat void the solution is appropriate for a crack with very sharp tip, 
i.e. a mathematical crack. The special case of longwave limit is given as 
an example. Since the solution given is analytic, computer display of 
results for different material systems and aspect ratios can easily be 
obtained, if needed. 
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